arXiv:1501.01634vl [astro-ph.CO] 7 Jan 2015 


Quintessence’s Last Stand? 


Eric V. Linder 

Berkeley Center for Cosmological Physics & Berkeley Lab, 

University of California, Berkeley, CA 94720, USA 
(Dated: January 9, 2015) 

Current cosmological data puts increasing pressure on models of dark energy in the freezing 
class, e.g. early dark energy or those with equation of state w substantially different from —1. We 
investigate to what extent data will distinguish the thawing class of quintessence from a cosmological 
constant. Since thawing dark energy deviates from w = — 1 only at late times, we find that deviations 
1 + w < 0.1 are difficult to see even with next generation measurements; however, modest redshift 
drift data can improve the sensitivity by a factor of two. Furthermore, technical naturalness prefers 
specific thawing models. 


I. INTRODUCTION 

Cosmic acceleration indicates crucial new physics ex¬ 
ists outside the standard model of particle physics and 
cosmology. Cosmological data imply the expansion his¬ 
tory of the universe acts like the matter component is 
supplemented with a fluid of strongly negative effective 
pressure: dark energy with an effective equation of state, 
or pressure to energy density, ratio w ~ —1. One possi¬ 
bility is the cosmological constant, with w = — 1 for all 
times, or redshifts z, and no spatial perturbations. For 
data sensitive to only the background cosmic expansion, 
or Hubble parameter H(z), the effects arc fully described 
by w(z ), even if there is no physical dark energy fluid but 
rather a modification of the form of the action. 

Cosmic microwave background (CMB) measurements 
impose two stringent constraints on the dark energy. 
From the distance to last scattering, the energy den¬ 
sity weighted time average of w(z) must be close to —1 
[1], From the structure of the CMB temperature power 
spectrum, the dark energy density around recombination 
(z ~ 10 3 ) is less than about 1% of the critical energy 
density [1-3] . Late time observations of the cosmic ex¬ 
pansion, such as supernovae or baryon acoustic oscilla¬ 
tion distances, also indicate that w ~ — 1 at recent times 
[4, 5]. These results disfavor the first quintessence mod¬ 
els, tracker models with insensitivity to initial conditions, 
but moreover most of the freezing region of dark energy 
dynamics (Ref. [6] showed that the effective dark energy 
field either began dominated by Hubble friction - thaw¬ 
ing fields - or by the steepness of the potential - freezing 
fields, unless fine tuned). 

The remaining dark energy behavior favored by data is 
the thawing class. While thawing fields do not have the 
trackers’ insensitivity to initial conditions, some of them 
do have the attractive quality of technical naturalness - 
protection against quantum radiative corrections. Since 
thawing models only deviate from w ~ —1 at late times, 
they are consistent with the data but are also harder 
to distinguish from a static cosmological constant, de¬ 
spite having very different physics implications. Given 
the current state of data, and the advancing plans for 
next generation experiments, we assess whether contin¬ 


ued consistency of the data with a cosmological constant 
A will enable all quintessence models to be disfavored, or 
at least relegated to a region fine tuned to be close to A. 

Beyond quintessence, cosmic acceleration models with 
deviations in growth of structure relative to their expan¬ 
sion behavior, such as from modified gravity theories, or 
clustered or coupled dark energy, can also be constrained 
by growth measurements. However within quintessence 
- canonical, minimally coupled dark energy - no new 
signatures arise within growth. Thus the question of 
whether we will be able to distinguish thawing dark en¬ 
ergy, and hence most of the remaining viable phase space 
of quintessence, from A will focus on the expansion his¬ 
tory. 

In Sec. II we discuss methods of treating the thaw¬ 
ing class of dark energy as a whole, drawing a clear dis¬ 
tinction with inappropriate slow roll assumptions. We 
examine the leverage of future data in separating thaw¬ 
ing quintessence from a cosmological constant in Sec. Ill, 
highlighting the role of potential redshift drift measure¬ 
ments. We explore some particle physics implications of 
fields nearly indistinguishable from A, and the virtues of 
technical naturalness, in Sec. IV and conclude in Sec. V. 

II. THAWING QUINTESSENCE 

The class of thawing quintessence includes common 
potentials such as monomials V((f>) ~ 4> n and pseudo- 
Nambu Goldstone boson (PNGB, or axion) fields with 
V{<j>) ~ [1 + cos(</>//)], where / is the symmetry break¬ 
ing energy scale. Two members of this class are of par¬ 
ticular interest since they have shift symmetry protect¬ 
ing against high energy radiative corrections. One is the 
PNGB case, long studied in terms of both natural in¬ 
flation [7] and dark energy [8]. The other is the linear 
potential [9, 10]. Recently the linear potential has gar¬ 
nered renewed interest in connection with both inflation 
[11] and modified gravity [12]. 

We are interested in constraining the thawing class as 
a whole, rather than individual models. To do this ef¬ 
ficiently it is useful to parametrize its behavior. The 
standard wo~w a form for the dark energy equation of 
state, w(a) = Wq + w a (l — a), where a = 1/(1 + z) is 
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the scale factor, works for a broad variety of dark energy 
models, beyond just the thawing class. Indeed it was de¬ 
rived from exact solutions of the field dynamics [13] (and 
is completely unrelated to a linear or Taylor expansion). 
The wq-Wo, form matches exact solutions of the observ¬ 
ables - distances and Hubble expansion H - to better 
than 0.1% for a broad range of currently allowed models 
[14]. However, because it involves two parameters, con¬ 
straints from even next generation data have difficulty 
distinguishing thawing models with wq < —0.9 from a 
cosmological constant at 95% confidence level, due to co- 
variances between parameters. 

Since we focus our interest on the thawing class, we 
could use more specialized parametrizations, such as [15- 
20]. We emphasize though that from an observational 
perspective the Wo~w a parametrization with its 0.1% 
matching is wholly sufficient. To obtain tighter model 
constraints, we would need a one parameter form. We 
cannot rely on the usual inflationary slow roll field ap¬ 
proximation since even thawing dark energy does not 
slow roll, in the sense that all terms in the Klein-Gordon 
equation for its dynamics are comparable [21]. For exam¬ 
ple, during the vast majority of e-folds when the field is 
frozen (e.g. during the matter dominated era), the terms 
are in the ratio of 1 : 2 : —3. Note also that a small field 
assumption is somewhat uncertain: while many thawing 
fields roll a distance A cf> « 0.24 Mpi to get to wo = —0.9 
today, they start at several Mpi from their minima. 

We can take two approaches to adopting a one param¬ 
eter thawing equation of state. Ref. [14] showed that 
not only was wo~w a an excellent approximation to ex¬ 
act solutions of the field evolution, and to observational 
quantities, but it could also be derived as a calibration 
relation for classes of dark energy dynamics. By appro¬ 
priate choice of wq and w a , the spread of evolution in 
the w-w' plane calibrated into narrow tracks for differ¬ 
ent dark energy physics. In particular, a broad range of 
thawing models could all be tightly fit by a constrained 
w a form: 

w a ~ -1.58 (1 + w 0 ) ■ (1) 

This combines all the virtues of the general wo~w a form 
(e.g. 0.1% distance reconstruction) with the leverage of 
a one parameter fit. 

While this fits the expansion observables (distances 
and Hubble parameter) superbly, an alternate approach 
is to follow the thawing physics more closely. For exam¬ 
ple, thawing models have w = — 1 at high redshift. This 
is not a problem for the form Eq. (1) since thawing dark 
energy density fades quickly into the past and so using 
w{z 1) = wo + w a ^ — 1 has negligible effect. (Recall 
that even accounting for this the distance to high red- 
shifts is accurately matched to better than 0.1%.) If we 
felt more comfortable approximating w(z) 1 despite it not 
being an observable, then the algebraic thawing form of 


Ref. [16] works quite well: 

/ i + b W-p / 3 

l+u,(a) = (l+u,oK , (2) 

where the constant b = 0.5. This is derived from the 
physics of how the field evolves upon leaving the matter 
dominated era and automatically has the correct high 
redshift behavior. 

We can turn it into a one parameter form by fixing 
p = 1, so the final algebraic form we will use is 

/ o \ 2 /3 

1 +w{a) = (1 +w 0 )a 3 (^ 1 + 2a3 J • ( 3 ) 

Note that of course it can fit a cosmological constant as 
well, with wq = —1. Both the (general) algebraic thawed 
and the general wo~w a forms have the advantage that 
the Hubble parameter H(z) is analytic. 

Figure 1 demonstrates the fit of the algebraic form to 
the exact solutions for thawing w(z ), with wq = —0.9. 
(Note that thawing models with greater present devia¬ 
tion from —1 are already in some tension with current 
data.) The agreement is better than 0.1% in w(z) (even 
better for the observables like distances and H ) for the 
favored models of the linear potential and PNGB, and 
better than 0.2% and 0.3% for the quadratic and quartic 
potentials. 



z 


FIG. 1. The one parameter algebraic thawing form (solid 
black curve, Eq. 3) accurately fits w(z) at all redshifts, for 
the thawing models of linear, quadratic, quartic, and PNGB 
potentials. 

While monomial potentials indeed have only one pa¬ 
rameter to describe the equation of state, PNGB models 
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have two: the steepness of the potential, given by the 
symmetry breaking scale /, and the initial field position. 
For steeper potentials, the algebraic form fit can degrade 
to « 1% in w (alternately we could adjust the values of b 
or p in the original algebraic form of Eq. 2, but there will 
still be a range of PNGB models that cannot be captured 
to much better than 1% in w(z) by only one parameter). 

Despite the success, it is important to note we are not 
trying to fit w(z) per se (as it is not an observable), but 
rather to obtain a realistic enough representation that 
our constraints on distinguishing thawing models from A 
in the next section are accurate. Table I demonstrates 
the success of both the algebraic thawer (Eq. 3) and con¬ 
strained wq-Wo, form (Eq. 1) in matching the observables. 
(Changing the value of, e.g., the matter density would 
slightly degrade the accuracy of fitting w(z), but improve 
the accuracy of fitting d or H since there is more freedom 
allowed for the fit.) The accuracy is more than sufficient 
for next generation observations seeking to distinguish 
thawing dark energy from a cosmological constant A. 


Model 

d(z) 

H{z) 

^lss 

Linear potential 

0.06% 

0.01% 

To 715 


0.08% 

0.02% 

10 -5 


0.16% 

0.03% 

10 -5 

PNGB {f/Mpi = 2) 

0.03% 

0.01% 

10 -5 

Constrained w a 

0.09% 

0.01% 

4 x 10“ 5 

A 

3.5% 

0.7% 

0.5% 


TABLE I. Maximum deviations in the observables - the co¬ 
moving distance d(z), Hubble parameter H(z), and distance 
to CMB last scattering di ss - over all redshifts, relative to the 
algebraic thawing form, are given for the exact solutions of 
various thawing dark energy models with wo = —0.9. The 
matter density and Hubble parameter are fixed (fi m = 0.3, 
h = 0.7). 


III. CONSTRAINING THAWERS 

As Table I shows, the maximum difference between 
thawing dark energy deviating from the cosmological con¬ 
stant at late times to wq = —0.9 is 3.5% in distance (at 
2 = 0.5). However, this was with fixing all other cos¬ 
mological parameters so in reality covariances make the 
distinction more challenging. Let us see if our two one- 
parameter thawing forms give data the leverage to dis¬ 
tinguish such dynamical dark energy from a cosmological 
constant. 

To test these models with expansion history measure¬ 
ments we consider supernova distances and the CMB dis¬ 
tance to last scattering. (We could use baryon acoustic 
oscillation distances, but supernovae are somewhat more 
sensitive to the equation of state, especially at the needed 
low redshift.) For supernovae (SN), we consider a future 
sample of 150 SN at z < 0.1, 900 between z = 0.1-1, and 
42 between z = 1 1.7, with a magnitude systematic of 
0.02(1 + z). Given the systematic control out to 2 = 1, 


this roughly represents a sample from ground based imag¬ 
ing surveys such as the Dark Energy Survey or Large 
Synoptic Survey Telescope (LSST). We also analyze an 
extended sample inspired by the DESIRE proposal [22] 
for the Euclid satellite, with 7200 SN between 2 = 0.1-1 
and 1000 at z = 1 1.6, with the same systematics con¬ 
trol out to z — 1.6. These are simply rough estimates 
of potential data. For the CMB we use Planck quality 
distance to last scattering measurement at the 0.2% level 
and a prior on the physical matter density D m /i 2 of 0.9%. 

The parameters for constraining the expansion history 
include the matter density D m , dynamical equation of 
state parameter wq, supernova absolute magnitude A4, 
and Hubble constant h. Spatial flatness is assumed. Us¬ 
ing the future data, the constraint on the equation of 
state parameter is er(u>o) ~ 0.06. If the thawing dark 
energy has reached wo = —0.9 today (and a stronger de¬ 
viation is already disfavored by present data), then this 
yields at most a ~ 1.7er distinction from A. A stronger 
discriminant would be useful. 

Interestingly, the extended distance sample does not 
improve substantially the constraints here. While the 
baseline sample gives tx(u>o) = 0.063 and 0.062 for the 
algebraic thawer and calibrated w a forms, the extended 
sample gives 0.056 and 0.055. The reason for this in¬ 
sensitivity is that thawers tend to deviate appreciably 
from A only at low redshift; for example thawers as rep¬ 
resented by the algebraic form have w{z = 1) = —0.98 
while eventually reaching w(z = 0) = —0.9. While su¬ 
pernova distances are indeed sensitive to the equation of 
state value at low redshift, they are not sensitive enough. 

The question then is what observational probe can 
more keenly test the low redshift equation of state. Re¬ 
cently, cosmic redshift drift has been recognized as highly 
sensitive at low redshift [23]. This constrains the quan¬ 
tity 


z = (1 + z) H 0 - H(z) . (4) 

For the (general) algebraic thawer the Hubble parameter 
can be written analytically as 



Ur 


-3 


+ (1 — D m ) e 


3(l + tr 0 ) 


[l-(aa 3 +/3) p/3 ] 


( 5 ) 

where a = 1/(1 + b), /3 = b/(l + b). Since redshift drift 
is wholly unproven, we conjecture a single, modest 5% 
measurement at redshift z = 0.5. This one data point 
improves the distinction between thawers and A by a 
factor of two, reducing cr(wo) to 0.030, marginalized over 
the other cosmological parameters, and hence a 3.3er dis¬ 
tinction from a cosmological constant. 


Of course as wo —» —1 the dark energy would be in¬ 
creasingly difficult to distinguish from A. Section IV dis¬ 
cusses particle physics implications of such near-A fields. 
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IV. FINE TUNED FIELDS 

Thawing fields have a limit in which they are identical 
to a cosmological constant, i.e. when wo = — 1, and so the 
full class can never be ruled out if the data are consistent 
with A. We examine whether fields in such a limit are 
particularly fine tuned or disfavored in particle physics 
terms. 

First consider quadratic and quartic potentials as 
thawers. Like all thawers, the fields do not roll very 
far during their evolution up to the present, as shown 
in Fig. 2. However, to attain even only wq = —0.9 to¬ 
day, the fields must start 2.6 Mpi and 5.1 Mp\ respec¬ 
tively from the minima. The potentials must be suffi¬ 
ciently steep locally to be able to reach ilde,o = 0.7 to¬ 
day once they are released from Hubble friction to roll. 
Values of u>o closer to —1 exacerbate the situation, with 
the fields needing to start at 7.9 Mp\ and 15.7Mpi re¬ 
spectively to reach, say, wo = —0.99 and hence be easily 
confused with A. Especially since such potentials have 
no protection against radiative corrections, we might be 
concerned about such superPlanckian values (despite the 
field rolling over less than a Planck mass). If we restrict 
the field initial position to lie within one Planck mass of 
the minimum, these potentials cannot deliver dark en¬ 
ergy density fide > 0.45 and 0.17 respectively. 



FIG. 2. The distance the field rolls up to the present is plotted 
vs the present equation of state for several thawer models. 
The field only traverses a short distance for, say, wo = —0.9 
but must start several Mpi from the potential minimum. 

The linear potential, however, is robust. For V(cf>) = 
Vo(l + a<^), while the slope of the potential aVo ~ 10 -120 
as for any quintessence, the parameter a itself is of order 
unity. For wq = —0.9 we have a = 0.72, and behavior 


close to A can be obtained without fine tuning a, e.g. 
wo = —0.99 comes from a = 0.25. Furthermore the lin¬ 
ear potential is substantially protected against quantum 
corrections. In the future it leads to a cosmic doomsday 
[10], actually essential for the recent solution of the cos¬ 
mological constant problem known as the sequester [12]. 
Short of doomsday, the linear potential can be extremely 
difficult to distinguish from A, however, and so upcoming 
observations will not be able to rule out fully the thawing 
class. 

For the PNGB model, a shift symmetry protects the 
potential form against quantum corrections. Neverthe¬ 
less, one might be concerned about the physics valid¬ 
ity if the symmetry energy scale / > Mp\. As / de¬ 
creases, and the potential steepens, it becomes harder 
to attain flde,o = 0.7. The field must be fine tuned to 
an initial position closer and closer to the potential 
maximum, out of all the range [0,7r/]. For fixed 4>i/ f , 
fide ,max oc / 2 , while for fixed fl de ,o, (<&//) max ^ ^ 
[14]. Thus small / is problematic. To obtain wq < —0.99 
with / = 0.5 Mpi requires faff < 0.18, i.e. only allow¬ 
ing 6% of the full range, while for / = 0.2 Mp\ a 0.5% 
fine tuning is required. Still, a reasonable region near 
f /Mpi ~ 1 is allowed and can give an expansion history 
effectively indistinguishable from A. 

Thus we see that nothing prevents thawing fields, and 
in particular the two technically natural and hence ro¬ 
bust models of the linear potential and PNGB, from ap¬ 
proaching A beyond the ability of upcoming experiments 
to distinguish. 


V. CONCLUSIONS 

Great strides have been made in constraining dark en¬ 
ergy physics with ever improving observations. Freezing 
dark energy, including early dark energy, is substantially 
limited by the data as an explanation for cosmic accel¬ 
eration. Probes of the expansion history will continue to 
become more incisive, and probes of the growth of large 
scale structure will impose tighter constraints on more 
elaborate theories involving dark energy clustering, cou¬ 
pling, and modified gravity. However even if the data 
continue to be consistent with a cosmological constant, 
quintessence remains a viable option in its thawing class. 

The standard dark energy equation of state 
parametrization w(a ) = wo + w a (l — a) is an ex¬ 
cellent global approximation over a wide range of 
freezer, thawed, or modified gravity models, to the 0.1% 
level in the observables. A constrained, one parameter 
form with w a = —1.58(1 + wq) can be used to model 
specifically the thawing class. Another option is to use 
the algebraic thawer form of Eq. (3) that follows the 
thawing physics to match w(z) as well as the observ¬ 
ables. We use these not as fitting forms per se, but to 
represent the thawing class as a whole for comparison to 
a cosmological constant, to evaluate the discriminatory 
leverage of forthcoming data. 
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We find that next generation distance data can reveal 
modest, but consistent distinctions of thawers from A if 
w 0 > —0.9. The leverage in distinguishing the physics 
doubles with inclusion of a single, 5% measurement by 
the prospective probe of cosmic redshift drift. This gives 
further motivation for its exploration and development. 
In particular, since the signature of the thawing deviation 
increases at very low redshift - e.g. the deviation is five 
times greater at z = 0 than at z = 1 then the recently 
identified power of redshift drift at low redshift [23] is 
especially valuable. 

While quintessence in the form of thawing fields can¬ 
not be ruled out even if the data become increasingly 
consistent with a cosmological constant, despite the very 
different physics, such behavior does point to preferred 
models. Both the linear potential and PNGB are tech¬ 
nically natural, easing high energy quantum effects by 


their shift symmetry. Indeed the linear potential has re¬ 
cently been found essential in one method for solving 
the original cosmological constant problem [12]. If fu¬ 
ture data show some definite signature in expansion or 
growth away from A, we will be able to focus our efforts 
on specific physical properties. Even if A continues to be 
a good fit, we will be drawn either to explaining its mag¬ 
nitude or to finding the informative high energy physics 
origins of these robust thawing quintessence potentials. 
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